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ABSTRACT 

We consider the effects of different criteria for determining where stars will form in gas on galactic 
scales, in simulations with high (1 pc) resolution, with explicitly resolved physics of GMC formation and 
destruction and stellar feedback from supernovae, radiation pressure, stellar winds, and photo-heating. 
We compare: (1) a self-gravity criterion (based on the local virial parameter and the assumption that 
self-gravitating gas collapses to high density in a single free-fall time), (2) a fixed density threshold, 
(3) a molecular-gas law, (4) a temperature threshold, (5) a requirement that the gas be Jeans -unstable, 
(6) a criteria that cooling times be shorter than dynamical times, and (7) a convergent- flow criterion. 
We consider all of these in both a MW-like and high-density (starburst or high-redshift) galaxy. With 
feedback present, all models produce identical integrated star formation rates (SFRs), in good agreement 
with the Kennicutt relation; without feedback all produce orders -of-magnitude excessive SFRs. This is 
totally dependent on feedback and independent of the SF law, even if the "local" collapse efficiency is 
100%. However, the predicted spatial and density distribution depend strongly on the SF criteria. Because 
cooling rates are generally fast within galaxy disks, and gas is turbulent, criteria (4)-(7) are very "weak" 
and spread the SF uniformly over most of the disk (down to densities n ~ 0.01 — 0.1 cm -3 ). A molecular 
criterion (3) localizes to slightly higher densities, but still a wide range; for metallicity near solar, it is 
almost identical to a fixed density threshold at n ~ 1cm -3 (well below the mean density in the central 
MW or starburst systems). A fixed density threshold (2) can always select the highest resolved densities, 
but must be adjusted both for simulation resolution and individual galaxy properties - the same threshold 
that works well in a MW-like simulation will select nearly all gas in a starburst. Binding criteria (1) 
tend to adaptively select the largest local over-densities, independent of galaxy model or resolution, and 
automatically predict clustered star formation. We argue that this SF model (possible with other secondary 
criteria) is most physically-motivated and presents significant numerical advantages in simulations with a 
large dynamic range. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — star formation: general — 
cosmology: theory 



1 INTRODUCTION 

Modeling star formation accurately is critical for any simulation of 
galaxy formation. However, cosmological and galaxy-scale simu- 
lations still cannot hope to resolve the spatial and density scales on 
which star formation actually occurs. As a result, simple "recipes" 
must be applied. For example, models typically impose some "lo- 
cal Schmidt law," where gas forms stars at a rate that scales as some 
power of the density; if the SFR per local free-fall time oc 1/ \[Gp 
were constant, this would be p* oc p/ta oc p 1 ' 5 , although in principle 
any other parameters can be used. 

Usually, applying these models alone would artificially spread 
star formation among all the gas in the simulation, even cosmo- 
logically pristine material at high temperatures. So some addi- 
tional criteria or restrictions must be included. Most commonly, 
this amounts to a simple density threshold: n > 0.1 cm -3 in many 
cosmological simulations. This is not to say such low-density mate- 
rial directly forms stars: rather, this corresponds crudely to densities 
where the thermal instability sets in, so some un-resolved fraction 
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of the material (which goes into the sub-grid scaling above) will be 
able to form stars. Other common requirements include restricting 
star formation to gas which is below some temperature, or Jeans 
unstable, or in convergent flows, or which has a short cooling time. 
Recently, various studies have considered molecular criteria: using 
some combination of density and metallicity to estimate a sub-grid 
molecular gas fraction and restricting star formation to the "molec- 
ular" gas ( [Robertson & Kravtsov|2"008l|Kuhlen et al.|2011| ). 

These criteria are not trivial. Star formation is observed to 
be highly clustered under essentially all conditions ( |Lada & Lada| 
2003 and references therein), and without some criteria such as 
those above, this is not captured in simulations. This does not just 
mean that stars are forming in the wrong places. "Smearing out" 
star formation over the disk dramatically suppresses the effects of 
stellar feedback ( Governato et al. 2010 ). Massive star clusters allow 
for e.g. overlapping SNe "bubbles" that can expand much more effi- 
ciently than individual SNe remnants. They also concentrate feed- 
back "where it is needed," i.e. it preferentially acts in the dense, 
star-forming gas. Strong radiation pressure effects arise when pho- 
tons are trapped in optical thick regions around embedded clusters 
(Hop kins et al.|201l) . And spreading out star formation leads to 
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spurious geometric cancellation between feedback sources. With- 
out a sufficiently strict minimum SF criterion, the ability of the gas 
to form realistic phase structure ( Saitoh et al. 2008), or blow winds 
that regulate it s baryon content ([Governato et al.|2010|) and form 
realistic disks (Pontzen & Governato 201 1 Gover nato et al.| 2004) 
can be fundamentally altered. 

Unfortunately, in practice the physical interpretation of these 
criteria often depends both on the resolved dynamic range of the 
simulation and on the mean properties of the galaxies being sim- 
ulated, and to obtain similar results they must be numerically ad- 
justed accordingly. 

However, it is increasingly clear that the ISM is governed by 
super- sonic turbulence over a wide range of scales. Consider, then, 
a locally self-gravitating region of the ISM "supported" by turbu- 
lence. In the absence of some feedback disrupting it or "pumping" 
the dispersion, the turbulent support will be damped in a single 
crossing time; as a result the region will collapse to arbitrarily high 
densities in about one free-fall time. This is true even if the en- 
ergy of contraction maintains a constant virial equilibrium at each 
radius (Hopkins 201 2bJ. At sufficiently high densities, eventually 
all of the above criteria must be met; so long as the SFR increases 
with density, eventually an order-unity fraction of the gas will be 
consumed into stars. So - in the absence of some self-regulation - 
the time-averaged SFR should be p* w p/f co iiapse ~ p/ts regardless 
of the "true" local star formation criteria/law. This is precisely what 
is seen in detailed simulations of turbulent media, in the absence of 
"pumping" to unbind collapsing regions CBallesteros-Paredes et al.| 
|2nTT||Padoan & Nordlund|20TT||Padoan et al.|2012) . 

We stress that this does not mean the total SFR, even within a 
dense parcel of gas, will actually be as large as pgas/fe. If feedback 
is present, it can self-regulate. As soon as some gas turns into stars, 
feedback can act and disrupt the bound material, terminating the 
star formation locally and suppressing nearby star formation even 
in the dense gas. In such a model, however, the "net" efficiency is 
actually predicted self-consistently from the feedback model, rather 
than imposed by the sub-grid model (by, say, forcing some by-hand 
low normalization of p'*(p). 

Of course, if a simulation resolves a bound region, then this 
collapse will be followed self-consistently. What we require is a 
mechanism to treat further collapse, where it would occur, below 
our resolution limit. In this paper, we propose a simple adaptive 
self-gravity criterion for star formation in galaxy-scale and cosmo- 
logical simulations, motivated by the numerical simulations above. 
We compare it to other common criteria in the literature, and ex- 
amine the implications for the equilibrium SFRs and both the pre- 
dicted spatial and density distributions of SF in different galaxy 
environments. 



2 A SIMPLE SELF-GRAVITY CRITERION 

On some scale Sr, self-gravity requires crl s + < f3 GM(< Sr) / Sr, 
where <r e ff includes the contributions from both rotational and ran- 
dom motions and f3 is an appropriate constant that depends on the 
internal structure in Sr. In practice the gas of interest is always 
highly super- sonic in our simulations (and in observations), so we 
can ignore the c s term here. We then obtain the usual virial param- 
eter 

a = a 2 eS Sr//3GM(<Sr). (1) 

To determine local binding, we wish to describe a in the limit 
where Sr is small. Then M(< Sr) — (47r/3) pSr 3 where p is the 



average density in Sr, and <r e ff — > Sr, or more formally 

a e 2 ft = p v (| V • v| 2 + | V x v| 2 ) Sr 2 = V (2) 

Here the V • v term accounts for the local radial velocity disper- 
sion and inflow/outflow motions, while the V x v term accounts 
for internal rotational/shear and tangential dispersion. The f3 v term 
depends on the internal structure again but is close to unity. 

Combining these terms, we can derive the formally resolution- 
independent criterion: 

Gp 

where /3' ~ 1/2 collects the order-unity terms above. This pre- 
f actor depends on the internal mass profile and velocity structure, 
but only weakly: for e..g a Plummer sphere or Hernquist (1990) 
mass distribution with pure isotropic, rotationally supported, or 
constant velocity gradient orbits the range is f3' w 0.5 — 0.6. 

This criterion is well-behaved in the local limit and does not 
explicitly depend on any numerical parameters of the simulation 
(spatial or mass resolution), nor does it require inserting any "ad 
hoc" threshold or normalization criterion. Because it depends only 
on the local velocity gradient and density, it is trivial to implement 
in either Lagrangian (SPH) or Eulerian (grid) codes (as compared 
to an explicit evaluation of the binding criterion over some resolved 
scale length, which requires a neighbor search, and can be pro- 
hibitively expensive in certain situations). 

Implicitly, the velocity gradients and average density are al- 
ways evaluated at the scale of the resolution limit. If going to higher 
resolution would change these quantities, then of course the cri- 
terion would give a different result (but as a consequence of the 
physical, not numerical, difference). 

3 THE SIMULATIONS 

The simulations used here are described in detail in Hopkins et al. 
( 201 1) (hereafter Paper I; see § 2 & Tables 1-3) and Hopkins et al. 
(2012b) (Paper II; § 2). We briefly summarize the most important 
properties here. The simulations were performed with the parallel 
TreeSPH code GADGET- 3 ( |Springel|2005] ). They include stars, dark 
matter, and gas, with cooling, shocks, star formation, and stellar 
feedback. 

3.1 Star Formation Criteria 

Star formation is allowed only in gas that meets some set of cri- 
teria, for example in density or temperature. Within the gas that is 
flagged as "star forming," our standard model assumes p* = ep/ts 
where ts is the free-fall time and e is some efficiency. Unless oth- 
erwise specified, we set e = 0.015, to match the average observed 
efficiency in dense gas (e.g.[Krumhorz & Tan 2007 and references 
therein). 

There are several criteria that can be imposed to determine the 
gas allowed to form stars: 

(1) Self-Gravity: We require a region be locally self- 
gravitating as described in §|2] i.e. a < 1. But these regions are 
assumed to collapse in a single free-fall time, so e = 1 (this approx- 
imates the results in individual cloud simulations of Padoan et al. 
|2012] >. 

(2) Density: Star formation is allowed above a simple density 
threshold n > no, where we adopt no = 100cm -3 to ensure this 
selects only over-dense gas inside of typical GMCs. 

(3) Molecular Gas: We calculate the molecular fraction fu 2 
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of all gas as a function of the local column density and metallic- 
ity following Krumholz & Gnedin (201 1) and allow star formation 
only from the molecular gas (i.e. multiply e by fn 2 ). 

(4) Temperature: We allow star formation only below a min- 
imum temperature T < Tmin- Here we adopt r m i n = \00K, (chosen 
to approach the minimum temperatures the simulation can resolve). 
At these temperatures, we expect this to be very similar to criterion 
(2). 

(5) Jeans Instability: We require the gas be locally 
Jeans-unstable below the resolution limit: c s < h sm \ \/4ttGp 
(h sm \ is the SPH smoothing length). Given the Lagrangian na- 
ture of the simulations, this translates to a temperature T < 
100K (mi/300 Mq)/ (/* sm i/10pc) where m,- is the particle mass. 

(6) Converging Flows: Star formation is allowed only in con- 
vergent flows, i.e. where V • v < 0. 

(7) Rapid Cooling: We allow star formation only in regions 
where the cooling time is less than the dynamical time, t coo \ < 

3.2 Cooling & Feedback 

Gas follows an atomic cooling curve with additional fine- structure 
cooling to rsj 10 K. At all the scales we resolve, the cooling time in 
dense gas tends to be much shorter than the dynamical time for any 
temperatures T > 10 4 K where the thermal pressure would be sig- 
nificant, and the minimum resolved scales are significantly larger 
than the sonic length. As a result, varying the cooling curve shape, 
magnitude, or metallicity dependence within an order of magnitude 
has no significant effect on any of our conclusions. 

Stellar feedback is included, from a variety of mechanisms. 

(1) Local Momentum-Driven Winds from Radiation Pressure, 
Supernovae, & Stellar Winds: Gas within a GMC (identified with 
an on-the-fly friends-of-friends algorithm) receives a direct mo- 
mentum flux from the stars in that cluster/clump. The momentum 
flux is P = PsNe + + Pmd, where the separate terms represent the 
direct momentum flux of SNe ejecta, stellar winds, and radiation 
pressure. The first two are directly tabulated for a single stellar pop- 
ulation as a function of age and metallicity Z and the flux is directed 
away from the stellar center. Because this is interior to clouds, the 
systems are always optically thick, so the latter is approximately 
Pmd ~ (1 + TiR)L inc ident/c, where 1 + tjr = 1 + S gas K m accounts 
for the absorption of the initial UV/optical flux and multiple scat- 
terings of the IR flux if the region is optically thick in the IR (with 
S gas calculated for each particle). 

(2) Supernova Shock-Heating: Gas shocked by supernovae 
can be heated to high temperatures. We tabulate the SNe Type-I and 
Type-II rates from |Mannucci et al.| ([2006 ) and STARBURST99, re- 
spectively, as a function of age and metallicity for all star particles 
and stochastically determine at each timestep if a SNe occurs. If so, 
the appropriate mechanical luminosity is injected as thermal energy 
in the gas within a smoothing length of the star particle. 

(3) Gas Recycling and Shock-Heating in Stellar Winds: Gas 
mass is returned to the ISM from stellar evolution, at a rate tab- 
ulated from SNe and stellar mass loss (integrated fraction ^0.3). 
The SNe heating is described above. Similarly, stellar winds are 
assumed to shock locally and inject the appropriate tabulated me- 
chanical luminosity L(t, Z) as a function of age and metallicity into 
the gas within a smoothing length. 

(4) Photo-Heating of HII Regions and Photo-Electric Heating: 
We also tabulate the rate of production of ionizing photons for each 
star particle; moving radially outwards from the star, we then ion- 
ize each neutral gas particle (using its density and state to determine 
the necessary photon number) until the photon budget is exhausted. 



Ionized sas is maintained at a minimum ~ 10 4 K until it falls out- 
side an HII region. Photo-electric heating is followed in a similar 
manner using the heating rates from Wol fire et al.| jl995). 

(5) Long-Range Radiation Pressure: Photons which escape the 
local GMC (not accounted for in (1)) can be absorbed at larger 
radii. Knowing the intrinsic SED of each star particle, we attenu- 
ate integrating the local gas density and gradients to convergence. 
The resulting "escaped" SED gives a flux that propagates to large 
distances, and can be treated in the same manner as the gravity tree 
to give the local net incident flux on a gas particle. The local ab- 
sorption is then calculated integrating over a frequency-dependent 
opacity that scales with metallicity, and the radiation pressure force 
is imparted. 

Details and numerical tests of these models are discussed in 
Paper II. All energy, mass, and momentum-injection rates are taken 
as-is from the stellar population models in STARBURST99, as- 
suming a Kroupa (2002 ) IMF, without any free parameters. Sub- 
tle variations in the implementation do not make significant differ- 
ences to our conclusions. Most important, we do not "turn off" or 
otherwise alter any of the cooling or hydrodynamics of the gas. 

3.3 Galaxy Models 

We implement the model in two distinct initial disk models, chosen 
to span a wide range in ISM densities. Each has a bulge, stellar and 
gaseous disk, halo, and central BH (although to isolate the role of 
stellar feedback, models for BH growth and feedback are disabled). 
At our standard resolution, each model has ^0.3 — 1 x 10 8 total 
particles, giving particle masses of ~ 500 Mq and typical ~ 5pc 
smoothing lengths in the dense gaQ and are run for a few orbital 
times each. The disk models include: 

(1) MW: a MW-like galaxy, with with baryonic mass Mb ar = 
7.1 x 1O 1O M (gas m g = 0.9 x 1O 1O M , bulge M b = 1.5 x 
10 10 Mq, the remainder in a stellar disk nid) and halo mass Mhaio = 
1.6 x 10 12 Mq. The gas (stellar) scale length is h g = 6.0kpc (hd = 
3.0). At standard resolution, gas particles and new stars have mass 
m « 5 00 Mq and the force softening e « 4pc. 

(2) HiZ/Starburst: a massive starburst disk with densities typi- 
cal of the central couple kpc in low-z galaxy mergers, or the larger- 
scale ISM in star-forming galaxies at z ~ 2 — 4. Here Mhaio = 1 .4 x 
10 12 Mq and baryonic (M bar , nib, ntd, m g ) = (10.7,0.7,3,7) x 
10 10 M with scale-lengths (hd, h 8 ) = (1.6, 3.2)kpc. 

4 RESULTS & DISCUSSION 

Figure[T]shows the total SFR versus time for the models, with feed- 
back enabled and disabled. We immediately see that the total SFR, 
hence the net SF efficiency, is almost totally independent of the SF 
criteria. However, it depends dramatically on feedback. With feed- 
back on, the global SFR is M* ~ 0.02Mg as /%n, in agreement with 
observations. Without feedback, it is M* ~ 0.5 M g as/% n [^] Regard- 
less of how mass turns into stars, a certain feedback strength, hence 
mass in young stars, is needed to inject enough momentum to offset 

Smoothing lengths are set adaptively as in |Springel & Hernq uist ( 2002), 
with an approximately constant 64 neighbors enclosed within the smoothing 
kernel. 

2 The apparent "convergence" of the feedback and no-feedback runs at late 
times occurs because the no-feedback runs exhaust all their gas in just a 
few dynamical times. Thus the SFR declines in absolute terms, even though 
it remains fixed in units of M gas /^dyn- This is shown explicitly in Paper 
I & Paper II. We also refer interested readers to Paper I, Fig. 1 1 for an 
explicit comparison of the simulations and the observed Schmidt-Kennicutt 
relation. 
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Figure 1. SFR vs. time for the simulations with different star formation 
prescriptions (as labeled), in the MW-like (top) and starburst/HiZ (bottom) 
disk models. The SF prescription has no effect on the actual SFRs (model 
differences, after the first couple dynamical times where non-equilibrium 
effects appear, are consistent with random variation). With realistic stellar 
feedback models, the SFR is entirely set by feedback (the mass in young 
stars needed to prevent runaway collapse). With no feedback, SF efficien- 
cies are extremely large (SFRs much larger than observed); with feedback, 
the SFRs agree well with observations in both models. In particular, despite 
having an instantaneous, local efficiency e = 1 in bound clouds, the actual 
average efficiency (e) of the "self-gravity" model is a realistic ~ 1 — 5%. 

dissipation and prevent runaway collapse. Because cooling is rapid, 
the galaxy can always find a way to get "enough" gas to the relevant 
densities/temperatures/binding criteria - at which point the SFR is 
a function of feedback efficiency alone. This is discussed in much 
greater detail in Paper I, where we show that, for example, chang- 
ing the local SF law (efficiencies, power-law indices, or density 
thresholds) makes no difference here, but changing the feedback 
strength directly changes the net efficiency. Without explicit feed- 
back models, a local SF model that turned self-gravitating gas into 
stars with an efficiency ~ 1 would grossly over-produce observed 
SFRs; however, once feedback is included, the mean SF efficiency 
even in dense gas is a self-regulating quantity and we can safely 
consider such a prescription^] 

The SF criteria does, however, affect the the spatial and den- 
sity distribution of star formation in each model, shown in Fig- 

3 In Paper II, we examine how different feedback mechanisms individually 
affect the SF efficiency. In the lower-density regime typical of the MW-like 
model, the dominant mechanism is a combination of photoionization heat- 
ing (and resulting warm gas pressure) and the momentum injection in over- 
lapping SNe explosions. In the higher-density regime of the HiZ model, it is 
predominantly radiation pressure on dusty gas from the nearly Eddington- 
limited starburst. 



MW Self-Gravity HiZ 
4 #. 




Figure 2. Gas surface density (intensity), with star-forming regions color- 
coded as red-yellow (with increasing specific SFR). The distribution of SF 
varies with SF law. Local binding criteria select locally over-dense regions 
in all cases. A fixed density threshold works well for most of the MW disk, 
but fails in the central, high-p regions of the HiZ model (where the mean 
density is above-threshold). A molecular law works reasonably well in the 
outer regions of the MW model, but "smears" SF among a much wider 
range of gas in spiral arms; in the HiZ/starburst nucleus model it iden- 
tifies all gas as molecular. Temperature, Jeans stability, cooling rate, and 
convergent- flow criteria select gas at nearly all densities. 
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Figure 3. SFR-weighted density distribution for the simulations with dif- 
ferent star formation prescriptions (labeled) in the MW-like (top) and star- 
burst/HiZ (bottom) disk models. The self-gravity criterion identifies the 
most dense regions. A fixed-density threshold simulation is dominated by 
SF near the threshold, which is much lower than the maximum densities in 
the HiZ model. A molecular criterion effectively corresponds to a much 
lower threshold; for the metallicities here « 1cm -3 . The other criteria 
spread the SF across almost all the gas, even at very low densities. 

ures |2|3| With arbitrarily high resolution, we would ideally want 
all the star formation to be highly spatially clustered, and concen- 
trated in absolute densities of > 10 6 cm -3 . Resolution limits mean 
this cannot be realized; the "next-best" aim of these star formation 
criteria is to identify the highest local over densities and highest 
resolvable densities, which will contain sub-regions that collapse 
further, and associate star formation with those regions. 

Our proposed self-gravity criterion (1) adaptively selects the 
over-dense regions in all galaxy models, giving a realistic "clumpy" 
and clustered morphology for the star formation. In the HiZ model 
in particular, it is able to associate star formation with a number of 
small sub-clumps in the central regions, with densities > 10 4 cm -3 . 
(Note that the apparent "bimodality" here is artificial, caused by a 
bottleneck at the highest densities allowed by our numerical reso- 
lution). 

A pure density criterion (2) works well when the threshold 
density is much larger than the background mean; however, it qual- 
itatively changes character when the threshold falls below the mean 
density (the central regions of the HiZ model). Then, obviously, the 
SF becomes "smeared" uniformly across all the gas, defeating the 
purpose of a density threshold in the first place. This is clear in 
the morphology in Fig. [2] where the "clumps" previously evident 
have been largely "smeared out" by more extended star formation. 



Likewise in Fig. [3] the SFR is dominated by gas near the thresh- 
old density, even though it is clear that densities up to 4 dex larger 
can be resolved. A number of authors have shown that this leads 
to unphysical distributions of stars and star formation, and can arti- 
ficially suppress the efficiency of feedback as well (see § [TJ. We 
stress that at least some of the cases here cannot simply be re- 
solved by raising the threshold further: for example, in simulating 
galaxy nuclei, the minimum density any clump must have simply 
to avoid tidal disruption (a clear requirement for star formation) 
scales oc r -3 , where r is the distance to the BH. Thus properly se- 
lecting "overdensities" at ~ 1 pc would require a threshold a factor 
of ~ 10 9 larger than the threshold at ^kpc - and the threshold we 
use on those scales is already much larger than the value chosen in 
most simulations! 

Identifying star formation with molecules (3) is reason- 
ably similar to the choice of a high threshold density or self- 
gravity, when we focus on regions of very low mean surface den- 
sity/opacity. When the average surface density is < IOMq pc -2 (at 
Z rsj Zq), the medium is not self- shielding and becomes atomic- 
dominated, so these criteria all similarly select overdense regions 
where rapid cooling has enabled collapse. But as soon as the den- 
sity rises much above this value, the (dense) gas is essentially all 
molecular, and the criterion becomes meaningless (distributing star 
formation equally among all gas). This "smears out" star forma- 
tion in the dense gas, evident in Figure [2] (which now appears as 
if it were effectively lower-resolution). In fact, for metallicities of 
> O.IZq, we see in Figure [5] that this criterion is nearly identical 
to invoking a relatively low "threshold" density of m ~ 1cm -3 . 
At lower metallicities significant differences appear (Kuhlei fet al.| 
|2011| , but there is almost no difference in an instantaneous sense 
between this model and a threshold density - the differences owe 
to the fact that at these metallicities cooling rates are sufficiently 
suppressed such that the cooling time is no longer short compared 
to the dynamical time ( [Glover & Clark|2012| >. 

A temperature criterion (4), Jeans criterion (5), and rapid cool- 
ing criterion (7) identify the "molecular" gas from criterion (3), but 
also include a wide range of gas at even lower densities. This can in- 
clude e.g. adiabatically cooled gas in winds and cold clumps which 
have been shredded by feedback. As a result SF in Figure|2]is effec- 
tively distributed over all the gas (clearly over several dex in Fig. [3] 
down to n ~ 0.01 — 0.1cm -3 ); the only reason it is concentrated 
towards the center at all is because of the global concentration of 
gas mass. 

An inflow/convergent flow criterion (6) has almost no effect 
on the distribution of star formation, except to randomly select ~ 
1/2 of the gas. This is because the medium is turbulent on all scales, 
so SIGN( V • v) is basically random. 

The binding energy criterion used here appears to have a num- 
ber of advantages over the traditional fixed density threshold crite- 
rion for following star formation in galaxy simulations. It is phys- 
ically well-motivated and removes the ambiguity associated with 
assigning a specific density - it also agrees well with the results of 
much higher-resolution simulations of star formation in individual 
local regions (Padoan et al. 2012 ). With realistic feedback models 
present, it is possible to set e = 1, i.e. have order-unity efficiency 
in bound regions, and correctly reproduce the observed SF efficien- 
cies both globally and in dense gas. This removes the dependence 
of the observed Schmidt-Kennicutt relation on all resolved scales 
on the "by hand" insertion of a specific efficiency. Interestingly, the 
total mass in dense (ft > 100 cm -3 ) gas is similar in the runs with a 
density threshold and e = 0.015, and those with a self-gravity cri- 
terion and e = 1. In the latter case, a large fraction of the dense 
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gas at any moment is not locally self-gravitating; the broad distri- 
bution of virial parameters in simulated GMCs is shown in Paper 
II, Figs. 17-18 & Dobbs e t al.|20li] So in this (local) sense as well 
as in the global average SFR shown here, the "average efficiencies" 
emerge similarly. A more quantitative calculation of the dense gas 
distributions as tracers of the local star formation efficiency is pre- 
sented in Hopkins et al. ( 2 012a| >. 

Most importantly, the self-gravity criterion is inherently adap- 
tive, and so allows for the simultaneous treatment of a wide dy- 
namic range - critical for simulations of e.g. galaxy mergers, ac- 
tive galactic nuclei, and high-resolution cosmological simulations. 
In these models, gas which is likely to be "star-forming" in one 
context (say a GMC in the outer parts of a galaxy disk) might be 
orders-of-magnitude below the densities needed for it to be even 
tidally bound in other regimes, so no single density threshold is 
practical for realistic simulation resolution limitations. For exam- 
ple, in the Milky Way, the central molecular zone is observed to 
have very high gas densities relative to the solar neighborhood and 
even many GMCs, but is not strongly self-gravitating, and so ap- 
pears to have a SFR far lower than what would be predicted by a 
simple density threshold argument (see Longmore et al.|2013] >. 

Of course it is possible to combine the criteria here, requir- 
ing self-gravity in addition to some density /temperature/molecular 
threshold, for example. However, based on our results, the addition 
of most of these criteria will not dramatically modify the results 
from a binding criterion alone. Moreover, it is not entirely obvi- 
ous if they add physical information - at low densities, for exam- 
ple, one might posit that a region should not form stars unless it is 
also cold. But if it does not form stars (and so cannot be disrupted 
by feedback), it would collapse (if it could be resolved) to much 
higher densities, at which point it should eventually become cold 
and molecular as well. The exception is the cooling time criterion 
- rapid collapse implicitly assumes efficient cooling; it is less clear 
what will happen to a clump that is bound but cannot dissipate. 
Likewise, one may wish to adopt other criteria for non-bound re- 
gions. For example, associating a low but non-zero efficiency even 
with un-bound but turbulent regions above some density thresh- 
old, to represent the fact that there is an (unresolved) distribu- 
tion of densities therein, some of which might be self-gravitating 
themselves (see e.g. |Krumholz & McKee|2005[ [Hopkins ||2012a| ). 
We have experimented with such a prescription (with e = 0.01 for 
ftcrit > 100cm -3 ), and find the contribution from the un-bound ma- 
terial is generally sub-dominant (but non-negligible). 
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Ultimately, testing which prescription is most accurate in sim- 
ulations should involve direct comparison with observations. Since 
we have shown that the different simulation criteria are degenerate 
in their predictions of the total SFR and star formation "efficiency," 
this is not a good observational constraint to test the models. How- 
ever, it is possible to measure how the observed SF in resolved 
galaxies is distributed with respect to the average gas density (or 
column density) on different resolved scales, essentially construct- 



ing a direct analog of Fig. p 
tial distribution of SF in Fig. 



The differences there and in the spa- 
. 2] should also be manifest in quantities 
such as the spatial correlation function of star formation (and young 
stars), and the sizes of different star-forming regions (e.g. the char- 
acteristic sizes of the star-forming portions of spiral arms or GMC 
complexes). The behavior of the local, small-scale star formation 
law (as a function of density or other parameters) also ultimately 
informs the "sub-grid" physics these models are intended to repro- 
duce. 
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